Impurity effect as a probe of the paring symmetry in BiS2 based superconductors 
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The impurity effects are studied for the BiS2- layered superconductors based on a two orbital 
model with the Bogoliubov-de-Gennes technique. The superconducting critical temperature (T c ) 
as a function of the impurity concentration is calculated. Significant reduction of T c is found for 
spin singlet nearest-neighboring paring d-wave state and the spin triplet next-nearest-neighboring 
(NNN) paring p-wave state, while no depression of T c for isotropic s-wave state. The single impurity 
effects for various paring states are also explored. The impurity resonance peak in the local density 
of states spectrum is found only for the d-wave state. For the spin triplet NNN paring p-wave 
state, two in-gap peaks occur in the case of positive impurity potential; and a single in-gap peak 
is found for the negative impurity potential, which shifts towards to lower energy with decreasing 
the potential strength. These results can be used to detect the paring symmetry of the BiS2- based 
superconductors. 
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Superconductors with layered crystal structures such 
as cupratesi, Sr 2 RuO^, MgB^, Na x Co 2 02^, and iron 
pnictides 5 - have generated enormous research interest. 
The low dimensionality can affect the electronic struc- 
ture and can realize high transition temperatures (T c ) 
and/or unconventional superconductivity mechanisms. 
Quite recently, superconductivity has been observed 
in the BiS2- based compounds 6 , such as B14O4S3 6 , 
LaOi_ x F x BiS 2 7 , NdOi_ x F x BiS 2 § , CeOi_ x F x BiSA and 
PrOo.5Fo.5BiS2 10 , where the BiS 2 - layer is a basic unit. 
Soon after that, considerable efforts have been paid 
on the investigation of the physical properties in this 
systemiir— . The discovery of the newly basic supercon- 
ducting (SC) layers may open new fields in physics and 
chemistry of low-dimensional superconductors. 

The BiS 2 - based materials compose of a stacking of 
BiS 2 - SC layers and the spacer layers, which resem- 
bles those of high-T c cuprates and the Fe-based super- 
conductors. From the first principles band structure 
calculations 6,11 , it predicts that the dominating bands 
for the electron conduction as well as for the super- 
conductivity are derived from the Bi 6p x and 6p y or- 
bits. According to the Hall effect and the magnetore- 
sistance measurements^, an exotic multi-band features 
are found in the B14O4S3 compounds and the SC pairing 
occurs in the one dimensional chains. Electrical resistiv- 
ity measurements under pressure^ reveal that B14O4S3 
and LaOi_ x F x BiS 2 have different T c versus pressure be- 
havior, and the Fermi surface of LaOi_ x F x BiS 2 may be 
located in the vicinity of some band edges leading to in- 
stability for superconductivity. In the CeOi_ x F x BiS 2 
system 9 , it is found that the parent phase is a bad 
metal, instead of a band insulator in the case of B14O4S3 
compounds 6 . By doping electrons into the system, super- 
conductivity appears together with an insulating normal 
stated, which is clearly different from the cuprates and 
the iron pnictides. 

The most fundamental issue is the superconductivity 



mechanism in this BiS 2 - based compound, which is still 
unclear. It is predicted that the electron-phonon coupling 
constant is as large as A = 0.85, which emphasizes that 
LaOi_ x F x BiS 2 is a strongly coupled electron-phonon 
superconductor—' 19 . While, a good nesting of the Fermi 
surface at wave vector k = (71", n, 0) has been foun d 11 ' 18 . 
It is also proposed that the SC pairing is strong and ex- 
ceeds the limit of the phonon mediated picture^. More- 
over, the anharmonic model calculation shows that the 
vicinity of the charge-density-wave (CDW) instability is 
essential for the superconductivity 18 , which is reminis- 
cent to iron pnictides whose proximity to spin-density- 
wave (SDW) is well established 5 . Thus, the spin fluctu- 
ation is also proposed to account for the SC pairing in 
this family^. 

Another puzzle is the SC paring symmetry of the BiS 2 - 
based superconductors. Based on the the first princi- 
ples calculation^, it is found that the Fermi level for 
the nominal composition is located in the vicinity of 
the topological change in the Fermi surface. This gives 
a possibility that the SC symmetry changes depending 
on the doping level, provided that the superconductiv- 
ity originates from the electric correlation. Considering 
various many-body interactions, possible paring states 
are proposed, such as sign conserving/reversing s-wave 
and d-wave state s 11 ' 17 . Furthermore, the first principles 
calculation indicates the quasi-one-dimensional bands 1 ^, 
which may leads to the spin-triplet pairing. In the 
CeOi_ x F x BiS 2 compounds, superconductivity and fer- 
romagnetism (FM) is found to coexist 9 , where FM may 
arise from the Ce moments. The coexistence of super- 
conductivity and FM challenges the spin-singlet paring 
mechanism, which reminds us of the Sr 2 Ru04 mate- 
rial where the short-range ferromagnetic spin fluctuations 
give rise to the triplet pairing with p-like symmetry^. 

The impurity effect is one of the most important 
tools for identifying the nature of the paring state and 
the microscopic properties, which has been success- 
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fully carried out in both conventional superconductors 27 
and unconventional ones, such as cuprates 2 ^ and iron 
pnictides22r— . In this paper, we will study theoreti- 
cally the impurity effect in the BiS2- based materials 
and look into the pairing symmetry based on the two- 
orbital model from the first principles calculation^ and 
the Bogoliubov-de-Gennes (BdG) equations. The impu- 
rity concentration dependence of the critical temperature 
is explored. Significant reduction of T c is found for spin 
singlet nearest-neighboring (NN) paring d-wave state and 
the spin triplet next-nearest-neighboring (NNN) paring 
p-wave state, while no depression of T c for isotropic s- 
wave state. The single impurity effects for various paring 
states are also calculated. The impurity resonance peak 
in the local density of states (LDOS) spectrum is found 
only for the d-wave state. The evolution of the reso- 
nance peak with the impurity strength is also calculated. 
For the spin triplet NNN paring p-wave state, two in-gap 
peaks occur in the case of positive impurity potential; 
and a single in-gap peak is found for the negative impu- 
rity potential, which shifts towards to lower energy with 
decreasing the potential strength. These results can be 
used to detect the paring symmetry of the BiS2- based 
superconductors. 

The starting model Hamiltonian with the hopping el- 
ements, the pairing terms and the impurity part is ex- 
pressed by 

H = H t + H A + H imp . (1) 

In the present work, we use the two orbital model with 
the hopping constants being from Ref. [11]. Thus, the 
hopping term H t can be expressed by 

H t = - E (*wV c 1 m <t c jV CT + H.c.) - <o E c U Ci ^> ( 2 ) 

where i,j are the site indices and /i, v = 1,2 are the 
orbital indices. And to is the chemical potential. Con- 
sidering the possibilities of both spin singlet and triplet 
pairings, the paring term is written as 

Here, ± is for spin-triplet and singlet pairings, respec- 
tively, and the pairing potential A^-^ is defined as 

A tjv = -%^(( c vT<>4) ± (Ci/4 c jvt))» where ViM" is the 
SC interaction. And -ffi mp is the impurity part of the 
Hamiltonian, written as 

ffimp = E V s c L^ c ^°- ( 4 ) 

i ra /X(T 

In this paper, both single- and multiple- impurity effects 
are studied, and only the intraorbital scattering by non- 
magnetic impurities is considered. 

Then, the Hamiltonian can be diagonalized by solving 



the BdG equations self-consistently, 

\ TA&, ) { v? vs v^ s ) > 

(5) 

where the Hamiltonian Hi ll j ua is expressed by, 

Hi^jua = —tiftjv — toSijS^v + V s 8i t i m . (6) 

The SC order parameter A,- w -„ and the local electron 
density (n,u) are obtained self-consistently: 

(7) 

(«*> = E Kt\ 2 f( E n) + E KA 1 f( E n)l (8) 
n n 

Here f(x) is the Fermi distribution function. The LDOS 
is expressed by 

Pi (w) = E[KVI 2 ^n - w) + K,\ 2 6(E n + u)], (9) 

where the delta function 8{x) has been approximated by 
r/-7r(a; 2 + T 2 ) with the quasiparticle damping T = 0.002. 

Following Ref. [17], we focus our studies on the optimal 
doped sample with x = 0.56 used throughout the present 
work. The paring strength is chosen the be Vi^jv = 0.6eV 
for all paring states. The numerical calculation is per- 
formed on a 24 x 24 lattice with the periodic boundary 
conditions. In this paper, the energy and length are mea- 
sured in units of eV and the Bi-Bi distance a, respectively. 
To calculate the LDOS, a 80 x 80 supercell technique is 
used. 

Based on our self-consistent calculation, various paring 
symmetries are obtained for the different paring states. 
For the spin singlet onsite paring, an isotropic s- wave 
state occurs. For the spin singlet NN paring, the paring 
symmetry is d- wave. To get the NN paring extended s- 
wave, we set Ai +X = Ai +y in the previous result. For the 
spin singlet NNN paring, due to the anisotropic hopping 
in the (x + y) and (— x + y) directions, the gap consists 
both s- and d- wave components, i.e., the s — d wave in 
the p x orbit and s + d wave in the p y orbit. To obtain 
the one dimensional spin triplet p- wave state, the paring 
interaction is only considered in the (x + y) direction of 
the p x orbit and in the (— x + y) direction of the p y orbit. 

We first study the multiple-impurity effect for various 
paring states. In this case, the impurities are randomly 
distributed in the 24 x 24 lattice and only the strong im- 
purity potential V s — lOeV is considered. The impurity 
concentration dependence of the reduced critical transi- 
tion temperature T c /T c q is presented in Fig. 1, where 
T c q is the critical transition temperature without any 
impurity. As seen, for the onsite paring s-wave state, 
the critical temperature T c does not decrease with the 
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FIG. 1: (Color online) The impurity concentration depen- 
dence of the reduced critical transition temperature T c /T c o 
for various paring states, where T c o is the critical transition 
temperature without any impurity. 



increasing impurity concentration. While, it decreases 
with increasing impurity concentration for the other par- 
ing states, such as spin singlet NN paring, NNN paring, 
and spin triplet NNN paring states. For the spin sin- 
glet NNN paring state, the critical temperature decreases 
rapidly in the low impurity concentration region (< 5%), 
then it almost saturates with further increasing the im- 
purity concentration. For the spin singlet NN paring d- 
wave state and the spin triplet NNN paring p-wave state, 
the critical transition temperature decreases with the in- 
creasing impurity concentration monotonically. While, 
the reduction of the critical temperature in the p-wave 
state is faster than that in the d-wave state. 

Generally, studying the nonmagnetic impurity effect 
on superconductors helps greatly to investigate the pair- 
ing symmetry. It is well known that superconductivity 
is robust against small concentrations of nonmagnetic 
impurities in conventional s-wave superconductors, ac- 
cording to the Anderson's theorem^. Since supercon- 
ductivity is due to the instability of the Fermi surface 
to pairing of time-reversed quasiparticle states, any per- 
turbation that does not lift the Kramers degeneracy of 
these states does not affect the mean-field SC transi- 
tion temperature. While, a magnetic impurity, owing 
to the effect of breaking the time reversal symmetry, 
can break Cooper pairs easily. For instance, MgB2 is a 
BCS-type supercondoctor with the electron-phonon cou- 
pling and its symmetry is a s-wave, revealed by enormous 
experiments^ 7 ., such as the isotope effect by the substitu- 
tion of n B2£ and the observation of a coherent peak in 
the nuclear magnetic resonance (NMR) experiment^ 9 .. Zn 
ions have a full d shell and are nominally nonmagnetic. 
Effect of Zn-substitution on the MgB2 superconductor 
shows that the critical transition temperature T c is al- 
most unchanged with small Zn-concentration, while T c is 
significantly suppressed by the substitution of magnetic 



ions^S. These results are consistent with the Anderson's 
theorem. For the BiS2— based superconductor, the den- 
sity functional based calculations shows that this mate- 
rial may be a conventional s-wave superconductor with 
large electron-phonon coupling constan t 18 ' 19 . If this is 
the case, the critical transition temperature T c should 
not be depressed by small concentration of non-magnetic 
impurities, as revealed by our numerical calculation pre- 
sented in Fig. 1. Note that there is an enhancement of 
T c by the nonmagnetic impurity substitution, which is 
perhaps due to the enhanced density of state near the 
Fermi level. 

In contrast to the isotropic s-wave superconduc- 
tors, nonmagnetic impurities are paring breakers in 
anisotropic superconductors, such as heavy fermion 
superconductors 41 , and cuprats^. Based on the sec- 
ond Born approximation and the strong-coupling Eliash- 
berg theory, impurity-induced T c suppression in d-wave 
cuprate superconductors has been predicted 4 ^. Ex- 
perimentally, Zn substitution for Cu in the high-T c 
cuprates dramatically suppresses superconductivity. For 
instance, it has been shown that an impurity concen- 
tration of 2-3 at.% (per Cu) reduces T c to half or less 
the value for unsubstituted systems over a wide region of 
hole doping level in Zn substituted La2- x Sr x Cu04 and 
YBa2Cu306+xr^— • In our self-consistent calculation, it 
is found that the paring symmetry is d-wave in the NN 
singlet paring states for BiS2— compounds. The critical 
transition temperature T c is suppressed significantly by 
small concentration of nonmagnetic impurities. This re- 
sult is similar to that of high T c cuprate superconductors, 
indicating the paring breaking effect in d-wave SC gap. 
While, in the NNN singlet paring state, T c is suppressed 
quickly for small impurity concentration then saturates 
with further increasing impurities. According to our cal- 
culation, the SC gap has both s- and d- wave compo- 
nents. Hence, the suppression of T c is the consequence 
of the paring break effect of the d-wave component in the 
SC gap. While, the saturation of T c is due to the s-wave 
component. 

In principle, nonmagnetic impurities are pair break- 
ers for higher-orbital momentum states, such as p-wave 
stated. The reduction of T c due to nonmagnetic im- 
purity scattering in Sr2Ru04 compound with p-wave 
like symmetry has been calculated from the conventional 
Abrikosov-Gorkov formula^, which has been proved by 
the substitution of nonmagnetic impurity Ti 4+ ion o 46 ' 47 . 
For the B1S2— based superconductors, possible spin 
triplet p-wave paring state is propose d 11 ' 17 . From our 
self-consistent calculation, the rapid depression of T c by 
nonmagnetic impurity substitution can be tested exper- 
imentally. 

We now turn to study the single- impurity effect in the 
LDOS spectra of various paring symmetries. The calcu- 
lated results are presented in Fig. 2(a)-2(e). Shown in 
Fig. 2(a) is the LDOS spectra around a nonmagnetic 
impurity for the isotropic s- wave paring state. As can 
be seen, the LDOS spectrum in the bulk has a U-shaped 
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FIG. 2: (Color online) LDOS spectra around a nonmagnetic 
impurity at various paring states. The impurity strength is 
set to be V s = lOeV. The solid (black) curves are the LDOS 
spectra in the bulk, the dotted (blue) curves are the LDOS 
spectra at the NN sites of the impurity, and the dashed (red) 
curves are the LDOS spectra at the NNN sites of the impurity, 
(a): Spin singlet of onsite paring state; (b): Spin singlet of 
NNN paring state; (c): Spin singlet of NN paring for d wave 
state; (d): Spin singlet of NN paring for extended s wave 
state; (e): Spin triplet of NNN paring for p wave state. 



bottom, indicating a fully SC gap. Outside the gap, the 
van Hove singularity peak is clearly seen. At the NN site 
of the impurity, the LDOS spectrum is almost unchanged 
compared with that in the bulk, while at the NNN site of 
the impurity, the SC coherence peaks in the LDOS spec- 
trum are depressed. From the above result, it is found 
that there is no impurity state in the isotropic s- wave 
paring state. In general, potential scattering impurities 
are not pair breaker in the s- wave case according to the 
Anderson theorem?-. Presented in Fig. 1(b) is the LDOS 
spectra for the spin singlet NNN paring state. The re- 
sults are similar to the case of the onsite paring state, 
i.e., the LDOS spectrum has a U-shaped bottom and no 
impurity state occurs inside the SC gap. 

The LDOS spectra for the spin singlet NN paring d- 
wave state are shown in Fig. 2(c). It is found that the 
LDOS spectrum has a V-shaped bottom in the bulk, in- 
dicating the nodal characteristics. At the NN site of the 
impurity, there is a sharp resonance peak at luq = 9meV 
in the LDOS spectrum. At the NNN site of the impu- 
rity, there is also a weak in-gap peak at ujq — — 9meV. 
Note that, the SC coherence peaks are also depressed ob- 



viously by the impurity at the NNN site. The resonance 
peak inside the SC gap at the NN site of the impurity is 
a general feature in a d-wave superconductor, which has 
already been reported in cuprate superconductors both 
theoretically^ and experimentally^. The LDOS spectra 
for the spin singlet NN paring s- wave state are presented 
in the Fig. 2(d). The LDOS spectrum in the bulk still 
has a V-shaped bottom. The LDOS spectrum at the NN 
site of the impurity is almost unchanged compared with 
that in the bulk. Similar to the spin singlet onsite and 
NNN paring states, there is no impurity state in this case. 
However, at the NNN site of the impurity, the SC coher- 
ence peaks are almost flattened, which is quite different 
from those in the previous cases. 

The LDOS spectra for the spin triplet NNN paring p- 
wave state are shown in Fig. 2(e). As one can see, the 
LDOS spectrum in the bulk is also V-shaped with nodal 
characteristics. Previously, a fully nodeless SC gap was 
found in the Sr2Ru04 superconductors from the scan- 
ning tunnelling microscopy experiments^, which is quite 
different with our result. This is most due to its complex 
paring symmetry with out nodes at the Fermi surface, 
i.e., smp x + isinpy wave or sin(p x + p y ) + ism(—p x +p y ) 
wave^S. At the NN site of the impurity, there is a small 
in-gap peak at uiq = 3mcV in the spectrum. At the 
NNN site of the impurity, two in-gap peaks are found at 
ujq = ±3meV, where the in-gap peak at the positive en- 
ergy is much stronger than that at the negative energy. 
Moreover, the SC coherence peaks are almost flattened 
by the impurity, similar to that of the spin singlet NN 
paring s-wave state. From the above discussion, one can 
concludes that there are impurity resonance peaks only in 
the spin singlet NN paring d- wave state. The resonance 
peak is very sharp in the LDOS spectrum at the NN site. 
While, two in-gap peaks occur in the spectrum at the 
NNN site for the p- wave paring state and the SC coher- 
ence peaks are almost flattened by the impurity. These 
results can be used to detect the paring symmetry in the 
BiS2- based superconductors. 

The potential strength response of the impurity states 
is also calculated. The impurity states at various poten- 
tial strength for the d-wave paring state are presented in 
Fig. 3(a)-3(e). As seen in Fig. 3(a)-3(b), for the pos- 
itive impurity potential the resonance peak at the NN 
site of the impurity decreases with decreasing potential 
strength. At the impurity potential V s — leV, the res- 
onance peak almost disappears and merges into the SC 
coherence peak. However, the situation is quite different 
in the case of the negative impurity potential. At strong 
impurity potential V s = — lOeV, there is also a sharp res- 
onance peak at the NN site of the impurity located at 
luq = 6meV, which is shifted to lower energy compared 
with that of the impurity potential V s — lOeV. With 
decreasing the impurity potential, as shown in Fig. 3(d), 
the resonance peak at the NN site of the impurity splits 
into two peaks located at ujq = ±4.5meV, and the peak 
at the negative energy is much stronger than that at the 
positive energy. At impurity potential V s = — leV, there 
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FIG. 3: (Color online) LDOS spectra around a nonmagnetic 
impurity at various potential strength for d-wave paring state, 
(a): V a = 3eV; (b): V s = leV; (c): V s = -lOeV; (d): V s = 
-2eV; (e):V s = -leV. 



model with the BdG technique. From the calculation 
of multiple- impurity effect, significant reduction of T c is 
found for spin singlet NN paring d-wave state and spin 
triplet NNN paring p-wave state, while no depression of 
T c for isotropic s-wave state. The single- impurity effects 
for various paring states are also calculated. The impu- 
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are also two in-gap peaks at the NN site of the impurity 
located at ujq = ±13.5meV, and the peak at the positive 
energy almost disappears. 

The potential strength response of the impurity states 
for the p-wave paring state is presented in Fig. 4(a)-4(c). 
For the positive impurity strength, the LDOS spectra are 
similar to that of V s = lOeV, i.e., two in-gap peaks oc- 
curring at the NNN site of the impurity. These results 
are not presented for simplicity. Here, we focus on the 
LDOS spectra of the negative impurity potential for the 
p-wave paring state. As can be seen, at the NN site of 
the impurity, the LDOS spectra are almost featureless 
compared with those in the bulk. However, there is an 
in-gap peak in the LDOS spectra at the NNN site of the 
impurity, and the SC coherence peaks are almost flat- 
tened. Moreover, the location of the in-gap peak shifts 
from u>o — 1.5meV at V s = — lOeV to ojo — — 4.5meV at 
V a = -1.5eV. 

In summary, the impurity effects are studied for the 
BiS2- layered superconductors based on a two orbital 



FIG. 4: (Color online) LDOS spectra around a nonmagnetic 
impurity at various potential strength for p-wave paring state, 
(a): V s = -lOeV; (b): V s = -3eV; (c): V s = -1.5eV. 



rity resonance peak in the LDOS spectrum is found only 
for the d-wave state. For the spin triplet NNN paring p- 
wave state, two in-gap peaks occur in the case of positive 
impurity potential; while a single in-gap peak is found 
for the negative impurity potential, which shifts towards 
to lower energy with decreasing the potential strength. 
These results can be used to detect the paring symmetry 
of the BiS2- based superconductors. 
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